14:4.1' 


'PE PRODUCT l OM 


NASA/TM- /f9 g 2q8282 


An X-ray pulsar with a 
superstrong magnetic field 
in the soft 7-ray repeater 
SGR1806-20 

C. Kouveliotou’", S. Dieters**, T. StrohmayerH, 

j. van Paradijsttl, G. J. Fishman', C. A. Meegan*, K. Hurley?, 

J. Kommer s=, I. Smith* , D. Frail** & T. Murakamitt 

• s ■ v, r Space Si u‘»ru.‘ /.rt/in/'iJiiiry, \.4.N4/A/aM/irt// Sjjftrc Fhgftl CrMfer. 

/ /;i *rf »»••//«* /5A‘/2. t/S.A 

1 ' 'n:«v» > fries Space Hr'ttinh A<>ocwtton 
Cnnct*uy of Alabama w Huntsville. Alabama J5.W, (AS’/t 
’ > l.abmaiory far High Tnergy Astrophysics. NA SA/Cn *dda rd Space Flight Center, 

I (nferthell. Maryland 207/1, USA 

; 1 Institute "/Uifmi Pnnnekock'' c- Center for Htgh •Energy 

A*:>op'iystCi. Kruttlatw 40 a fOOti 5/ Amsterdam, Ihc Netherlands 

1 1 Sptuv Si inters {.aborntnry, University of Californio. Rcrkclcv. 

■ Cf/tf’nriff 9-/7.?<> 7450, U.S/l 
» 

j * Massachusetts Institute of Technology. MS .'7-62?. Cambridge, 

Masmehusetn 02130. USA 

i Rice University, Space Physics and Astronomy, 6100 South Main. MS 108 , 
i Hon nun. Texas 77251-1892, (75.A 

| “ STAG. VIA. 100 } f.opczvillc Rond. Socorro. New Mexico 3780 / , C/S/. 

1 • ** ISAS. 2-1-1 Yoshiitnriai, Sngamihnnr, Kano gown, japan 229 

j 

Soft v-ray repeaters (SGRs) emit multiple, brief { —0 . 1 -s) intense 
outbursts of low-energy 7 -rays. They are extremely rare 1 — 
j 1 href 3 * 4 are known in our Galaxy and one* in the Large Magellanic 
- Cloud. Two SGRs are associated 5 "' with young supernova rem- 
i mints (SNRs), and therefore most probably with neutron stars, 

‘ but it remains a puzzle why SGRs are so different from "normal’ 

; radio pulsars. Here we report the discovery of pulsations in the 
i persistent X-rny flux of SGR1806 - 20, with a period of 7.47 s and 
: a spindown rate of 2.6 x 10 " 3 s yr * 1 . Wc argue that the spindown 
I is due to magnetic dipole emission and find that the pulsar age and 
I (dipolar) magnetic field strength are -1,500 years and 
j 8 x I0 M gauss, respectively. Our observations demonstrate the 
: existence of ‘magnetars’, neutron stars with magnetic fields about 
! 100 limes stronger than those of radio pulsars, and support earlier 
j suggestions*'* that SGR bursts are caused by neutron-star ‘crust- 
j quakes' produced by magnetic stresses. The ‘magnetar’ birth rate 
[ is about one per millennium — a substantial fraction of that of 
i 

j Figure 1 Fast-Fcu^r-transform power spectrum computed lor the combined 
1 RX TF /PCa observations o f 5/6. 7 . 13 and 1 7 November 1936 in me 2-24 keV range. 

! data on ie November contained too m*ny thirsts to be useful) The no iviau a 1 
| powers have oeen normalized by the average noise evei to puny the distribution 
I oi cuse (>cwer? mro aporoximate agreement with The x 3 distribution with 2 
| degrees oi treeoom Tne oata were sampled m 0 r>-s bins The 0 i3375-Hz pulsed 
1 s*yn;ii a°n ts tvst na-monic a; 0 26760 Hr siano out wen apove the noise level 
j TrotJaia naveoecn b.vycentre cc^ectcd {usng the co-ordinates of the LOv star 
[ mentio-eo m the i*!xt 'or the source position! 
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radio pulsars. Thus our results may explain uln muhc vMTs Im vi- 
no radio pulsars. 

SGR180<S - 20 became extremely active between Oct , 1 Wi 
and November 1997, when over 40 intense bursts and numerous ' 
weaker ones were defected' 0 with the Burst And Trans ion: Source ! 
Experiment (BATSE) on board the Compton Gamnu-K.iv Ohser- 1 
vatory (CGRO). V\V observed SCR 1806 - 20 with :1 k* X-Rav 

Timing Explorer (1LXTE) rive times between s and 18 Non ember , 
1996. starting five days after the first triggered hutsi deletion with 
BATSE. (Information on the archival data from RX7F ?c.A and 
ASCA is available dt http://heasarc.gsfc. na>a gov ) Dtr.r.e these 
observations", the source emitted .series of outbursts in .1 b..nchmg 
mode, never seen before. The intensity of the outbursts .w well ns 
the ‘bunching' mode, varied significantly: mitu-outburs:- were 
interlaced with very intense ones and the rate of bur>o varied 
from bunch to hunch (5. Dieters e/ «/., manuscript m piopj- iuon). 

We made a period search of rhr data after excluding ah hursrs 
from the time series. The data were then energy-selected for 2- < 
24 keV X-rays, background subtracted and binned at (1.5 > resolu- 
tion. The resulting light curve was searched for periodicities 
between 0.03 and I Hz. by calculating a fast-Fourier transform 
power spectrum (Fig. 1 ). The peaks in the spectrum are centred 
on the fundamental frequency of 0. 13375 H/ (period of ~.4”655s) ; 
and its first harmonic at 0.26750 Hz. We find no significant power in 
any other frequency in the searched range. The probabihtv that we I 
detect a signal at the fundamental frequency this strong b\ chance 
coincidence is I x io 1 (taking into account the number of trials, 

1.9 X I0\ and the probability per trial, 5X10“ %0 ). 

To determine the fundamental period, alt data sets were then 
corrected to the Solar System baryccntre and separately folded at the 
longest detected period of 7.47655s, and sub harmonics thereof. 
These sub-harmonic folds showed multiple 7.47655-s pulses, which 
were identical, within statistical errors. We therefore conclude that 
7.47655 s is the fundamental period. 

To refine our period estimate, wc have combined all Promotional 
Counter Array (PCA) data and used a boot-strapping* method, where 
the best period obtained from the shortest data span is used :o count 
cycles and estimate the period over the next longest span ofda:^. .uul so 
forth until a consistent period is found over the entire data -e: This is 
achieved bv cross-correlating the phase folds from data sublet > u>in« 
our best period estimate and fitting the resulting cross-correlation < 
peaks with a gaussian. The phase offset is then converted tv* a time 
difference, a new period is calculated and the procedure N con- 
tinued for the next (longest) span of data. Our errors in the period 
and period derivative were estimated by performing Monte Carlo 
simulations of 25 data sets, whereby wc added gaussian noise with 
the data variance to the bins of each phase fold. 



N A J U IlL | VOi 793 I 2 I MAY I . 


235 


\ 4 - -14 


mhture ppglo: : 10 ^ 


P.03 


I 1- ,v 1hV-1 


letters to nature 

V\V fine! that j contain period docs i wn hr the pulse arrival rimes 
well . a iuI we derived an ephemeris of the pulse arrival times from the 
minimum m the phase-folded pulse light curve (sec Fig. 2): 

r„„„(TjD) = r„ + np +'~ pw' 

! where 7A is the reference rime in truncated Julian days 
j (TJD ~ ID - 24*40000.5), P is the pulsar period, P is the period 
dern alive and ;V is the number of elapsed cycles since T 0 . We find 
thj! T. = 10392. 45007L037(I)(TFD), P = 7.47655 1(3) s and 
P - 2.84 ! .4) x 10 '• ss * 

1 Using this ephemeris, the corresponding pulse profile of all the 
’ UX'f ii/EC.A data is shown in Fig. 2a. The full amplitude of the pulsed 
i tiux is 1.0 r 0 2 countss ' (between 2 and 24 keV). To estimate the 
i point source Count rate, we have folded the spectrum derived* 2 from 
I the observations of the Japanese satellite ASCA of the persistent X- 
ray source identified'' with SGR1806 - 20 through the detector 
response of the RXTE/PCA. It predicts 3.9 counts s" 1 for a point- 
like source detected with PCA in the 2-24 keV range; the corre- 
sponding pulsed flux amplitude is — 26%. 

I We have searched the archival ASCA data of SGR 1806 - 20 for 
i the presence of a periodic signal over a period range that encom- 
passed the original period and possible extrapolations given by our 
P measurement. A significant signal was found in the 1993 data set 
(epoch 3 JD 9275.85) at 7.4685 1 ± 0.00025 s with full amplitude of 
i pulsed flux “~23%. The probability of this peak to appear by 
| chance in the ASCA data is 3.6 X 10 4 , (taking into account the 
; number of trials. 377, and the probability per trial, 9.66 X 10 ~ 7 ). 

| Figure 2b shows the pulse profile of the ASCA folded light curve. We 
1 have a marginal detection of the pulsation in the 1995 data set 
| U'.473S * 0.00 1 s, epoch TJD 10006.72). The ASCA and RXTE/ 
i PC.A periods arc consistent with a long-term average spindown 
j rate. P = (S3 - 0.3) x W lt ss‘ l . The ASCA detection of the 
! periodicity clearly associates the period with SGR1806 - 20 rather 
I than with a hitherto unknown pulsar within the 1° field of view of 
j rhe RXTE/PCA. 

The period we have found in SGR1806 - 20 is very similar to the 
i S.0-s period found 3 in SGR0526 — 66 during a three-minute 
’ interval following its extremely strong 5 March 1979 burst. This 
| strongly suggests that the mechanism that produces SGR events is 
1 associated with a particular type of slowly rotaring neutron stars. 
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Figure 2 L;;och folded pulse proves of SGR 1806-20 a RXTE/PCA daw in 20 
onase o.ns a: the 'undamentai oer.od ol 7 4/655; s. Oorween 2 and 24keV. b. 
ASCA aota (from ’he ’993 observation of SGR 1 806 - 20) m 12 phase bins (to 
increase ‘he s«gnoi to no-se ^atio) at ire fundamental period of 73685:25s 


Whar is rhe energy source of the SCR iiuiewcmv ami bur>r 
emission? The (unabsorbed) llux U-lOkcV) iH t> 1c quiescent 
emission from SGR 1806 — 20 measured 17 with ASCA is 
1 X 1 0 11 ergem “ : s ’ ’. Assuming isotropic emission and a 
distance" for SCR1S06 - 20 of -14 kpc, the 2-10 keV luminosity 
of the source is « 2 x 10" ergs ‘ 1 . Similarly, the 

SGR 1 806 - 20 burst peak luminosity" (25-60 keV) corresponds 
to -10 41 ergs' 1 , which is ^ 10’ times higher than the Eddington 
luminosity of a neutron star of 1.4 solar masses (1.4 A/ Two Kpes 
of energy sources have been proposed to explain the emission from 
SGR s: accretion, and energy released by the decay of rhe magnetic 
field of a strongly ( ^ 1 0 1 J G ) magnetized neutron star, a ^-called* 
‘magnetar'. (Rorational energy loss, £ roi . is not sufficient to explain 
the observed X-ray luminosity of SGR 1806 - 20 For the observ ed P 
and P we find - 10" ergs' 1 two orders of magnitude below 
^sck) As we argue below, the 'magnetar’ model is the onlv one that 
can account for the observed properties of SGR 1806 - 20. 

SGR1806 - 20 has been identified with the SNR G 1 0 0 - 0.3, 
whose radio morphology is that of a plcrion'* (a synchrotron nebula 
powered by the relativistic wind of a young pulsar ); rhis suggests 
that it is centrally powered by a compact source. The X-rav 
spectrum 13 of the ASCA persistent source is a power law with 
photon index 2.2, and hydrogen column densitv, 
N h = 6 X 10 33 cm ~ 3 . It does nor show the emission lines com- 
monly seen in SNRs, which shows that the contribution from 
shocked gas is relatively small and further supports the idea that 
it is an (X-ray) plerion. In addition, radio observations" of 
SGR 1 806 - 20 show varying spatial structure ar an angular scale 
of several arcseconds, perhaps a jet. All the above observations 
indicate that the SNR is powered by a wind of relativistic particles 
Thompson and Duncan* estimate that the particle luminosity, L 
from SGR 1 806 - 20 is of the order of 1 0' 7 ergs" 1 . 

The coincidence of SGR1806 - 20 with an LBV (ref. 18). a rare 
type of evolved massive star, raises the possibility that the SGR 
persistent X-ray luminosity is caused by accretion from the LBV 
stellar wind. Accretion-powered emission, however, can be excluded 
by the following argument*, for reasonable (L. Kaper, personal 
communication) parameters of the LBV (wind velocity 
500 kms“\ mass loss rate ~ 10~‘ M 0 yr - \ radius -100 
solar radii, and the distance a between the two stars 2* l0 L 'cm): 
at the accretion radius (r 4Ci — 2GM NS VJ : , where is the neutron • 


b 



Phase 

Between 2 and 10 keV The plots in the two panels cover iwc pe'-ods io r c' 3 '*y 
The ASCA Dulse profile is nor derived by extrapolation o' he oxtc/pca per»od 
and period der.vative. it is ?ound by independent ocnod searches I 
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csMitv o.cr u d lw rhc t cl.it iv i>nc wind power: ng the 
'V,., “ /-,,/f‘iTTi is 1.2 x m dyncm \ This 

pic»ore exceeds the siclkir wind r.im pressure (P, 9n , “ 

; iJ,. v r ;. V' 25 dyn cm where p w is the density in the wind) by 
■ .itmoNt two orders of magnitude. Thus, even if the SGR source 
were m binary svMem with the 1.11V, accretion is highly unlikely to 
< »ccur 

I he radio ddui, liowevcr, truv provide evidence for the orbital 
• period of the system, that could be confirmed with a long term 
monitoring of rhc source with RXTE. Frail et at." report that the jet* 
like structure associated with the plcrion has rotated by 50* over an 
: interval of 1.4 yr. It is possible that this structure is the result of a 
magnetotail ot SGR 1806 — 20 produced by the strong wind of the 
LBV, I n such a case, wc would expect the system orbital period to be 
| of the order of 1 0 vr. 

Tins leaves us with the magnetar model. Because in this model 
j rotational losses do operate on the neutron star, we can estimate 
[ the spindown pulsar age, r, and magnetic field, B, from its 
i period and period derivative, according to 19 : r = Pf 2p and 
1 n - 3.2 x 10 respectively. Substituting the above values 
i for P and for the long-term P we find r-»l,500yr and 
j b =* 8 x I0 M G. Thompson and Blaes 20 have shown that these 
’ formulae arc of approximate validity only when the particle lumin- 
osity ftom a neutron star is much lower than its standard inagnetic 
j dipole luminosity. In particular, when the particle pressure is 
mgh enough to exceed the dipole pressure inside the light cylinder 
of the neutron star, the magnetic field lines will he combed out 
beyond a radius R y , where R v /R^ - (fliKi $ tV2£ p ) ,, \ Here B n and 
are the dipole magnetic held and the neutron star radius, 
respectively. This raises (he magnetic field strength at the light 
1 cylinder, and therefore the magnetic dipole torque. Following their 
: analysis, wc find for the magnetic field and spindown age of 
SGR1806 20, —2 X 10 N G and — 8,000yr, respectively. The 

: totter value is in good agreement with an earlier age estimate 1 * of 
— 10.000 yr for G !0.0 *— 0.3 based on angular diameter versus 
; surface brightness argument. 

The magnetic field value is the highest measured so far (its 
uncorrected value is similar to that estimated 2 ' for GB790305 and 

‘ for :: if- 184) - 045 in the young SNR Kes 73); it supports the model 

| for SGRs developed by Thompson and Duncan". According to these 
j workers, the SGR bursts are triggered by cracking of the neutron- 
star crust caused by magnetic stresses, which leads to sudden 
injection ot Alfv£n waves into the magnetosphere, particle accel- 
j eraiion. and formation of an optically thick pair plasma. Very high 
: magnetic fields would explain the super-Eddington luminosities of 
I SGRs, hv suppressing the electron scattering cross-section of 
j P ,u, f<>ns by a factor *£T’. In addition, the decay of the magnetic 
I field heats the neutron star interior, giving rise to persistent thermal 
j soft X-ray emission from the surface. Superstrong fields also lead to 
rapid spindown of a neutron star, which after a time interval of the 
I order of 10 yr will rotate at a period in the 10-s range. Thus, unless 
G 10.0 -0.3 is not n pierionic SNR, the conclusion that 
• SGR 1806 - 20 is a magnetar seems unavoidable. 

; SGRs have several striking simiki rines with a small group of X-ray 
I pulsars (only -6 sources are known so far), called ‘anomalous’ X- 
; iay pulsars (AXPs)*- ; \ which arc very different from normal binary 
i X-ray pulsars, with respect to their spin period distribution (all 
| between 6 and 12 s). their very soft X-ray spectra, and the absence 
j of evidence for a binary companion. Several AXPs are associated 
’ wuh young SNRs; their very small distances to the galactic plane 
i shows that AXPs are likewise recently formed neutron sta is' 4 . Non- 
detection of orbital Doppler shifts in their X-ray pulsations indi- 
cates that if they have companions, these must be very low-mass 
‘ stars. 


magnetic 
'■uli time r 


that the «ige range oi AXPs > Invited, .uni . ,n 
dipole spindown the nj.uh ptiunf nines rather Nio 
as f' : . For a particular model of ihe tf-fieM de^av. Thompson anti 
Duncan^ hnd for the period l> ^ rf . : (rf 1(>' yr) and 

the X-ray luminosity L v =*= 5 > I b'Tr/IO 4 yr) ‘ : ‘ erg s * (rhc latter 
expression depends somewhat on details nr the neutron star 
equation of state) Here U,. - B/IO'-'G, R„ = K v> /lir C m. arid 
M,. 4 =M ns / 1.4 M. The observed X-rav luminosities and spin periods 
therefore follow naturally from the model. 

According to a previous 1 analysis, the total number ut SGRs m 
out galaxy is < 7. If the age estimate of SCR 1806 - 2n t s tvpical 
for these objects, their birth rate is of the order of one per 
millennium, that is. -10% of the birth rate of normal radio 
pulsars* 07 . Neutron stars with superstrong magnetic fields may, 
therefore, constitute a non-negligrble fraction of new|v formed 
neutron stars and account for at least part of the discrepancy 
between the birth rates of core-collapse supernovac and of radio 
pulsars''*. ft is unclear whether these magnetars represent the tail of a 
very broad neutron-star magnetic -field distribution, or whether 
they form a separate class of young neutron stars with >upcr$rrong 
magnetic fields. 

We point out that the value of ptP of SGR I 806 - 20 is smaller 
than any of the values (only five measured so far) for AXPs. This 
may suggest that SGRs arc the earliest phase in c he evolution of 
magnetars, followed by a phase in which they appear as AXPs. □ 
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j The accretion model of AXPs fails to explain the narrow dis- 
j trihution of their observed pulse periods and X-ray luminosities. In 
tlie magnetar model the observed limited period range reflects 
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